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The theory for condensation of higher fermionic clusters is developed. Fully selfconsistent nonlin- 
ear equations for the quartet order parameter in strongly coupled fermionic systems are established 
and solved. The breakdown of the quasiparticle picture is pointed out. Derivation of numerically 
tractable approximation is described. The momentum projected factorisation ansatz of Ref. |2l[ for 
the order parameter is employed again. As a definite example the condensation of a particles in 
nuclear matter is worked out. 



I. INTRODUCTION 

Quartet condensation is relevant in several domains of 
physics. It is so far mostly considered in nuclear physics 
with the strongly bound alpha particle cluster playing 
a dominant role in certain states of lighter nuclei [l[ 
and, eventually, also in the surface of heavy nuclei as 
may be indicated by observed alpha decay processes. Al- 
pha particle condensation may eventually occur in com- 
pact stars m. However, trapping of multi component 
fermionic atoms makes fast progress. The case of trions 
is already quite advanced, experimentally 0, Q and also 
with theory [Ml. One may trap fermions with four dif- 
ferent 'colors' in the near future. This is a prerequisite 
for quartet formation and quartet condensation. Theo- 
retical work on this subject already has appeared [l3 - [l8j . 
The condensation of bi- excitons in semi-conductors also 
may be of relevance [l^, [H . 

In the recent past we already have published several 
papers on alpha particle condensation in infinite nuclear 
matter. For instance, we studied the onset of alpha parti- 
cle condensation and evaluated the corresponding critical 
temperature [2l|, [2^ with a procedure analogous to the 
pairing case via a four body generalization of the famous 
Thouless criterion However, quartet condensation 

not only is formally much more complicated than con- 
densation of (Cooper) pairs, it turns out that also certain 
aspects of the physics are quite different. The most strik- 
ing feature being that quartets only exist in the so-called 
BEC (Bose-Einstcin condensate) limit where they do not 
overlap very much in space. Contrary to the pairing case 
quartets cannot strongly intermingle in real space and, 
therefore, a coherence length much longer than the inter 
alpha particle distance cannot exist. We will give the 
reason for this different behavior in the main text. 

In this paper we will treat quartet condensation at 
zero temperature and establish and solve a full nonlin- 
ear equation for the quartet order parameter which will 
be the analog to the gap equation for pairing. A di- 



rect solution of such a highly non linear four body prob- 
lem seems hopeless, even in homogeneous matter. How- 
ever, we recently showed that a very simplifying approxi- 
mation works very well, at least in the domain of neg- 
ative chemical potentials, i.e. in the strong coupling 
regime. This approximation consists in making a mean 
field ansatz, i.e. a Slater determinant, for the quartet 
but projected onto zero total momentum as it is relevant 
for condensation plj . Such a mean field treatment may 
work each time the quartet is in its lowest energy config- 
uration. However, one may think of generalisations for 
excited configurations as well. Since the four fermions 
of the quartet are all different in the nuclear case (pro- 
ton/neutron, spin up/down), all are in the same OS mean 
field wave function and the problem boils down in the 
end to solve the equation in iterating on this single one 
particle mean field wave function. The same happens, of 
course, if the four fermions are components of an F = 3/2 
spin as in [16]. The problem still is complicated but can 
be solved with effective interactions of separable form. 
It was shown that that approximation gives comparable 
results as with the four body Faddeev-Yakubovsky equa- 
tion with a more realistic interaction in the nuclear case. 
So, in the present work we calculate the order parame- 
ter equation with the projected mean field ansatz and a 
separable potential. 

The paper is organized as follows: in the next section, 
we discuss the BCS gap equation. In Sec. Illli we show 
the full expression of the single particle mass operator 
for quartet condensation. Since this is too difficult to 
calculate numerically, we suggest an approximate mass 
operator in Sec. IIVI Before we show the numerical result 
employed by the approximated mass operator, we discuss 
in Sec. |V]the significant difference between pairing con- 
densation and quartet condensation through the different 
level densities involved in the condensation processes. In 
Sec. Ell we present the results. Finally we conclude in the 
last section. In the Appendix, we describe the detailed 
derivation of the equations and discuss various methods 
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to formulate the order parameter equation for quartet 
condensation. 



II. RECAPITULATION OF THE PAIRING 

CASE 




In order to prepare the terrain for our procedure in 
the quartetting case, in this section we want to rederive 
standard BCS theory in a way somewhat differing from 
the usual. 

The one body Green's function (GF) for BCS is repre- 
sented by p6| 



Sll' 



where Mf.f^{uj) is the BCS mass operator 



(1) 
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where 



^i2 = -iE 



1'12,34(C4C3) 



(2) 
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with (cc) being the thermal average of the pair operator. 

The indices 1,2,3,... correspond to momentum and 
spin. In nuclear matter, we have to add isospin. W12.34 
is the antisymmetrized matrix element of the two-body 
interaction (wi2,34 = -^12,43 = -W2i,34)- The single 
particle energies Si are in principle given by the kinetic 
energies plus the mean field shifts. The direct term is 
in homogeneous matter a constant which can be incor- 
porated into the chemical potential and the Fock-term 
gives rise to an effective mass. Since we will mostly deal 
with very low density nuclear matter, we will not con- 
sider a mass renormalization here. Therefore, in Eq. ([!]) 
we have ei = fc^/(2m) — fii with the chemical potential 
fii which, in principle, contains the direct part of the 
mean field. We have attached an index on the chemical 
potential, since in principal it can depend on the various 
fermionic components involved. However, in this work 
we always will consider fully symmetric situations and 
henceforth we will suppose that the chemical potentials 
of all fermionic species are equal and, therefore, drop the 
index. 

Fig. [T] is the graphical representation of the BCS mass 
operator. As shown in Eq. ([2]) and Figure [TJ the BCS 
mass operator consists of the 2 particle-1 hole (2plh) 
GF between two-body vertices factorized into two order 
parameters and the free Ih GF. 

On the other hand, the in-medium Schrodinger equa- 
tion for the order parameter is of the following form 



(C2C1) = - 



1 



Pi 



£1 +£2 



XI ^"12,1'2'(C2'C1'), (4) 



FIG. 1: Graphic representation of the BGS mass operator in 
Eq. © 




FIG. 2: Some diagram for the mass operator. 



where pi is the occupation number derived from 



Pi 



du , 
2^ 



2ImGi;i(cj + i77)/(w) 



(5) 



with the Fermi distribution function /(w) = [e"/-^H- 1]"-'^ 
and a positive infinitesimal of ry, as indicated. 

In the standard BCS theory, pairs in time reversed 
states are considered, i.e. in Eqs. ^ and (U) taking 2 = 
1, p\ = pj, and El = Ej, where 1 is the time reversal 
conjugate of quantum numbers 1. For Eq. ([T]), we obtain 
the imaginary part of the one body Green function as 



— IniGi;i(a; + 177) = 



(6) 



with Aj^j = Ai and E\ ~ -y/e^ + A^. This is equivalent 
to solving the usual gap equation at finite temperature as 
can easily be deduced from the spectral function obtained 
from © [illii]: 



V 



(7) 



Note that Eq. ^ resembles a particle-particle RPA 
equation (26j with, however, renormalized occupation 
numbers. One could, therefore, also consider Eq. (jH) as a 
single pole approximation to the so-called renormalized 
RPA, well known in the literature, see, e.g., Ref. [27| . 



III. SINGLE PARTICLE MASS OPERATOR 
AND QUARTET CONDENSATION 

Formally it is straightforward to generalize the pair- 
ing case to the quartet case. A typical diagram of the 
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mass operator needed for quartet condensation is shown 
in Fig. [21 As seen, to make the quartet order parameter 
appear in the single particle mass operator, we need to 
express it by the 4p3h GF. In order to get to it, we must 
successively open phase space from one particle to 2plh, 
to 3p2h, to 4p3h. This needs three interaction vertices 
on each side as shown in Fig. [21 Details of the derivation 
will be given in the Appendix [Bl We obtain for the Ip 
mass operator with quartet condensation 



rm4;5678(^)(c8C7C6C5) 



72/3/4 + f2hf4 

to + £234 
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where summation convention over repeated indices is to 
be understood, and 



^123; 
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The effective four body vertex F^''^ , evaluated perturba- 
tivcly to third order in the interaction is 



^ 1234;5678V"^^ 

r(3) (, ,\ 

^ 12'3':4'5'6' \^ ^ 
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2 W - £3'4' + £2' 



^W3'2,45^2'3<54'6 



+ 2"4'3,45'^2'2^3'6 ~ 2 ^2'6,23'53'4<54'5 



(11) 



where £i23 -- = £i + £2 + £3 H , and /j = 1 - fi with 

the Fermi distribution function fi ~ fi^i)- 

One may wonder why in (jlOp and (jlip no Fermi func- 
tion factors appear with the propagators. This, however, 
is a general feature of vertices coupling lower configura- 
tion spaces to higher ones. A fameous example in the 
literature is given by coupling Iplh and 2p2h spaces as 
this appears in the damping of zero sound modes, see, e.g. 
[29| . One may notice that the absence of Pauli blocking 
factors opens up phase space and, therefore, enhances the 
coupling. 



One can ask the question whether in Eq. ([8]) the un- 
correlated mean field 3h GF should be used. One may 
think to include the hole GF's into the selfconsistent cycle 
or even include higher correlations. One should notice, 
however, that in BCS this is not done and in Eq. ([2|) the 
mean field Ih GF is used. This has a good reason, since 
BCS theory is based on a variational wave function which 
fully respects the Pauli principle. Should a self consistent 
hole GF be used in ([2]), this property would be lost. We, 
therefore, also stick to the mean field hole GF's in the 
quartet case. 

We also need the in-medium four-body Schrodinger 
equation for the order parameter, in analogy with the 
pairing case shown in Eq. Q in the previous section. It 
is given by 



£l234(c4C3C2Ci) 

+ ^ Vi234;l'2'3'4'(c4'C3'C2'Ci') 



(12) 



1'2'3'4' 



where 



1234;1'2'3'4' 



= {I- Pi- P2)-jVi2-V2'533.'5ii' 
+ (1 - Pi - P3)^Wi3;i/3'(522"544' 



(13) 



-|- permutations. 

Details are given in [22j and Appendix [Cl 

We consider symmetric (nuclear) matter. In this case, 
we can give a fully symmetric order parameter of ex- 
change between two particles with respect to momenta: 

(C4C3C2C1) ^ h,MM,k.^o^ (^4) 

where the spin-isospin singlet wave function is repre- 
sented by and we consider here a spin-isospin inde- 
pendent two-body interaction: 



«12,34 -> %fc,,fc3fc,v 



j('5siS3^tlt3^S2S4'5t2t4 
— 5siSi5tit4.5s2S35t2t3) (15) 

with Si (ti) of spin (isospin) index. Vj:^!:^ ^ is symmet- 
ric with respect to exchange of the momenta: v^^^^ ^^^^ = 



^felfe2,fe4fe3 ■ 



Then, Eq. (|T^ is explicitly written as 



E 



n 



\ (277)3 

(l-p(fcl)-p(fc2))ffc^fc,,fc,fc, 

x(27r)3,5(fc3 - fc;,) (277)3,5(^4 

+ (1-P(fcl)-P(fc3))%fc3,fe^fe^ 

x{2TTfd{k2 - k'2){2TTf6{k4 
-^permutations 



(16) 



A sketch of the quartet mass operator is shown in Fig. [31 
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FIG. 3: Graphic representation of the mass operator for alpha 
condensation in Eq. ((8)1. 




FIG. 4: The graphical representation of the approximate a- 
BEC mass operator M"!""'^* of Eq. 



IV. APPROXIMATE QUARTET MASS 
OPERATOR 



From Eqs. (jlOp and ([TT|) it becomes evident that the 
fuh evaluation of the Ip f-^ 4p3h vertices are too com- 
plicated to be evaluated exactly. However, with quite 
reasonable approximations, one arrives at a numerically 
manageable expression. Since this discussion involves 
lengthy, quite technical details, we relegate it to Ap- 
pendix |D] and only give the final result here: 

^quartet j-^-j ^ ^ Ai234(/2/3/4 + /2/3/4)A^234 ^rj-^ 



In Eq. (jl8|) . we put a renormalization factor A' in front 
of the vertex which shall effectively account for the ap- 
proximations considered, see App.|Dj It is, however, very 
fortunate that the final result is independent of the value 
of A' and, therefore, one also may drop it. This some- 
what surprising effect is due to selfconsistency and self- 
readjustment of the solution. It is demonstrated in the 
Appendix |E] for the simple case of ordinary pairing. 

The complexity of the calculation still is further much 
reduced for the order parameter (|14p with our mean field 
ansatz projected on zero total momentum, as already 
very successfully employed in [2l[ 



h^Mk = ¥'(|fcl|)V'(|fc2|)v(|fc3|)¥'(|fc4|) 

x{2TTfS{h+k2 + k3 + k4). (19) 
It should be pointed out that this product ansatz with 
four identical OS* single particle wave functions is typi- 
cal for a ground state configuration of the quartet. Ex- 
cited configurations with wave functions of higher nodal 
structures may eventually be envisaged for other physi- 
cal situations. We also would like to mention that the 
momentum conserving 6 function induces strong correla- 
tions among the four particles and is, therefore, a rather 
non trivial variational wave function. 

For the two-body interaction of wj^g^ g^jt^ in Eq. (|15p . 
we employ the same separable form as done already in our 
previous publication on the quartet critical temperature 
in Ref. (ill: 



234 



OJ + £234 



fcl - fc2 , , ks - k4 

>'W[ )w( ) 

X {2n)6^{ki +k2-k3- (20) 



where the quartet 'gap' matrix A1234 is given by 

A1234 = A'-Ui2,l'2"533"544'(ci'C2'C3'C4'). (18) 

The graphical representation of the approximate M'^"'''*"* 
is shown in Fig. 3] 



with the form factor w{k) = w{\k\) = e"*"'^/^^. An exam- 
ple for strength and range parameters are given in Sec. lVIl 
below. 

With these simplifications, the mass operator ([TT]) is 
independent of spin and isospin, and therefore it can be 
reduced to the following four-fold integral 



Mi"^''"='(fci,a;) 



1 



(4^^)^ Jo 



dKK^ I dti I dkk' 



' J{\K ~ fciD/df + fc|)/(|f - k\) + f{\K fciD/df + fc|)/(|f ~ fc|) 



X («^(|fci-f l))'(^(lf + fc|))'(¥'(lf -fc|))' 



dk'k'^ 

r 



.K 



.K 



dt'w{k')^{\- + k'\M\--k'\) 



(21) 



with ti ^ {K-ki)/{Kki), t = {K-k)/{Kk), and t' ^ {K ■ 



k') / [Kk'). Here we represented the Fermi distribution 
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function as fi ^ /(^i) = /(^i)- In this expression any 
strength factor of the vertices has been dropped, see our 
above argument and Appendix |e1 For the imaginary part 
of this expression an energy conserving delta function 
comes instead of the full denominator and then the 4D 
integral can be reduced to a 3D one. How this goes in 
detail is again explained in Appendix |F] The real part 
of quartet then obtained from the imaginary part via 
a dispersion integral: 

ReM^"^'"='(fc,w + i?7) 



i-oo W - 



where V denotes the Cauchy principal value. 



The occupation numbers are finally obtained from 



p{k) = - J ^2ImG(fc,a; + *r,)/H 

du} -2ImMi"^''*°'(fc,w + m) 

= — fi^)- (23) 

27r [uj - e{k) - RcA/q"''''*<'*(fc,a; + i?/))^ + (ImMq"'i''*«t(fc, w + iri)Y 

I 

The equation for the order parameter (10) is formally tion. For completeness we again give the equations for 
not changed from Eqs. (4)- (7) of [2l| but the occupation the single particle wave function ip{k) 
numbers are calculated selfconsistently with above equa- 



A{k)ip{k)+W{k) + ?,C{k)Lp{k)=Q, (24) 

where 



{2„Y 



v'(k')v'(P'), (25) 



x^{k[)ip{k2Mk'^)^^ik3)^^{k^){2nf5{k + k2 + h + h) 
= J dKK^ j ^dt{\-p{k)-p{P))w{Q)^{P) J dk'k'' J ^dt'wiQ'Mk'MP') 

X J dk"k"^ dt"^^{k")^^{P"), (26) 

f d^k2 d^k3 d^ki d^k'n d^k'o ,^ , ,k2-ks^ rk2-k'n ^ 

m = j J^J^J^J^J^i^~Pik.)-piks))Xwi^)wi^)i2nm^^ 

Xip{k2)ip{k'^)ip{k3)'p{k'^)ip^{k^){2TTf5{k + k2 + h + h) 



J dKK' J dt^^P) j dk'k'' J^^ dt'{l p{k') p{P'))w{Q'Mk'MP') 
xjdk"k"^j dt"w{Q")ip{k")ip{P"). (27) 



3J 

1 r rl 



with 



P = 



VK^ + fc2 + 2Kkt, 



P' = y^K'^ + fc'2 + 2Kk't\ 
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P" 


= \Jk'^ + 


fc"2 + 


2Kk"t", 


Q 


= ^K^/A 


+ fc2 - 


\- Kkt, 


Q' 


= ^JK^I^ 


+ fc'2 


+ Kk't', 


Q" 


= VK^i 


+ fc"2 


+ Kk"t" 



As mentioned, in these equations the occupation 
numbers p{k) shah be calculated selfconsistently from 
Eq. 

Because of its particular importance, before the presen- 
tation of the results, we first will discuss in the following 
section the 3h level density. 



V. THREE HOLE LEVEL DENSITY 

In what follows a crucial role will be played by the 
3 hole propagator entering the mass operator. Since its 
influence on quartet condensation will be radically differ- 
ent from the corresponding Ih propagator in the pairing 
case, we will pay special attention to it in this section. 
In mean field approximation we can write 



iIlMl±llMA = G^''^\k,MM,^), 



£l23 



(28) 



where at r.h.s. we dropped spin and isospin indices since 
we consider unpolarised (nuclear) matter. We immedi- 
ately see a strong difference with the pairing case. There, 
only a single hole line enters whose numerator is (see pre- 
vious section) fi + fi = 1, see Eq. ^ Therefore, no single 
particle occupation numbers appear in the numerator of 
a single hole propagator. This difference between the 
three hole and one hole propagators leads to strong con- 
sequences. This is best demonstrated in Fig. [5] with the 
three hole level density 



d^ki d^k2 d^ks 
(27r)3 (27r)3 {2nf 
d^ki d^k2 d^ks 



lijiG^^''\ki,k2,k3;uj + iT]) 



where 



(27r)3 (27r)3 (27r)3 
x(/l/2/3 + ./l/2./3)7r(5(w + ei + £2 + £3) 

^^jl^^[dKK^ 
x(f{k)F{K,p) + f{k)F{K,p) 

kmax = v^2m(3^ - u), 



v/4m(3/i-a;) - 2fc2, 



3myLt — mu! 



k^ 
T 
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(29) 



(30) 
(31) 

(32) 




-40 -20 20 40 60 80 100 
uj (MeV) 



FIG. 5: 3h evel densities defined in Eq. (|29|) for various values 
of the chemical potential at a zero temperature. 



and 



F{K,p) = 
FiK,p) = 



dtf{\^-^+p\)fi\^-m, (33) 
rfi/(lf +Pl)/(lf -Pi) (34) 



with t = K -p/i^Kp). 

In Fig. Owe show the level density at zero temperature 
(/(w) = 6'(— w)), where it is calculated with the proton 
mass m = 938.27MeV (natural unit). Two cases have to 
be considered, chemical potential /i positive or negative. 
In the latter case we have binding of the quartet. Let 
us first discuss the case /i > 0. We remark that in this 
case, the 3h level density goes through zero at lo ~ 0, i.e. 
just in the region where the quartet correlations should 
appear. This is a strong difference with the pairing case 
where the Ih level density does not feel any influence 
from the medium and, therefore, the corresponding level 
density varies (neglecting the mean field for the sake of 
the argument) like in free space with the square root of 
energy. In particular, this means that the level density 
is finite at the Fermi level. This is a dramatic difference 
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FIG. 6: Single particle wave function ip{k) in fc-space (left), for 
r-space ip{r) (middle), and occupation numbers (right) at jj. — 
—5.26 (top), —1.63 (middle) and 0.55 (bottom). The r-space 
wave function <^(r) is derived from the Fourier transform of 

ip{k) by ip{r) = / d^fce** >(fc)/(27r)^ The dashed line in the 
left panels correspond to the Gaussian with same norm and 
rms momentum as f{k). 
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FIG. 7: -ImM''""-''*'=*(fci, w + ir?) in Eq. dm as a function of 
for = -4.9MeV (left) and for fj, = 0.55MeV (right) at zero 
temperature. 



with the quartet case and explains why Cooper pairs can 
strongly overlap whereas for quartets this is impossible 
as we will see below. We also would like to point out 
that the 3h level density is just the mirror to the 3p level 
density which has been discussed in psl ]. 

For the case where < there is nothing very special, 
besides the fact that it only is non- vanishing for negative 
values of UJ and that the upper boundary is given by 
UJ ~ 3/i. Therefore, the level density of Eq. is zero 
for UJ > 3fi. 



VI. RESULTS AND DISCUSSION 

At first let us mention that in this pilot application of 
our selfconsistent quartet theory, we only will consider 
the zero temperature case. As a definite physical exam- 
ple, we will treat the case of nuclear physics with the 
particularly strongly bound quartet, the a particle. It 
should be pointed out, however, that if scaled appropri- 
ately all energies and lengths can be transformed to other 
physical systems. For the nuclear case it is convenient to 
measure energies in in Fermi energies ep = 35MeV and 
lengths in inverse Fermi momentum kp^ ~ 1.35~^fm. 

To determine the order parameter, the calculation it- 
erates in the following cycle: 
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FIG. 8: -2ImG(fc,a; + i-q) in Eq. ^ as function of u) for 
fi = -4.9MeV (top) and for ^ = 0.55MeV (bottom) at a zero 
temperature. 



1. Fix a chemical potential ^ 
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2. Give an initial trial order parameter ip{k) 

3. Calculate the mass operator Af'5"^'^*°*(fc, w) of 
Eq. dH]) with if{k) 

4. Calculate the occupation numbers p(k) with 

^quartet^-^^^) from Eq. 



5. Substituting ip{k) and p{k) into Eqs. pS]). p6)) . and 
([?7| , derive the "new" single particle wave function 
c/3(fc) from 



-3B{k) 



Aik)+^C{k)' 



(35) 



6. Quit the cycle once ip{k) has converged, while go 
to 3. otherwise. 

The single particle wave functions and occupation 
numbers obtained from the above cycle are shown in 
Fig. m We take A = -992MeV fm^ and b = 1.43fm-i to 
get the binding energy of the free a-particle (— 28.3MeV) 
and its radius (1.7fm). We also insert the Gaussian wave 
function with same rms momentum as the single particle 
wave function in the left figures in Fig. [51 As shown in 
Fig. |6l the single particle wave function is sharper than a 
Gaussian. There is the term A{k) of Eq. in the de- 
nominator of (p{k) in Eq. (pSj) . and A{k) involves the fac- 
tor (fc^/(2m) — 4/i). Hence, <z'(fc) is closer to a Lorentzian 
rather than to a Gaussian 

We could not obtain a convergent wave function for 
H > 0.55MeV. This difficulty is of the same origin as in 
the case of our calculation of the critical temperature for 
a particle condensation. This stems from the fact that 
for larger positive values of the chemical potential, the 
denominator of ^{k) in Eq. psp at a certain value of k 
becomes zero, while the numerator is finite. In case of 
pair condensation, as shown in Eq. (|4]), both the denomi- 
nator and numerator become zero at the same value of k, 
this being a further crucial difference between the quar- 
tet and pairing cases. In the r.h.s. panels of Fig. [S] we 
also show the corresponding occupation numbers. We 
see that they are very small. However, they increase for 
increasing values of the chemical potential. For fj, = 0.55 
MeV the maximum of the occupation still only attains 
0.35 what is far away from the saturation value of one. 
What really happens for larger values of the chemical 
potential, is unclear. Surely, as discussed in Sec. |V] the 
situation for the quartet case is completely different from 
the standard pairing case. This is due to the fact, as al- 
ready mentioned, that the 3h level density goes through 
zero at = 0, i.e. just at the place where the quar- 
tet correlation should build up for positive values of fi. 
Due to this fact, the inhibition to go into the positive 
fi regime is here even stronger than in the case of the 
critical temperature [2l| . 

The situation in the quartet case is also in so far much 
different, as the the 3h GF produces a considerable imag- 
inary part of the mass operator. Figure[7]shows the imag- 
inary part of the approximate quartet mass operator of 



Eq. ([2T|) for /X < and p > 0. These large values of the 
damping rate imply a strong violation of the quasiparti- 
cle picture. In Fig. [5] we show the spectral function of the 
single particle GF. Contrary to the pairing case with its 
sharp quasiparticle pole, we here only find a very broad 
distribution, implying that the quasiparticle picture is 
completely destroyed. How to formulate a theory which 
goes continuously from the quartet case into the pair- 
ing case, is an open question. One solution could be to 
start right from the beginning with an in medium four 
body equation which contains a superfiuid phase. When 
the quartet phase disappears, the superfiuid phase may 
remain. Such investigations shall be done in the future. 



VII. CONCLUSION 

We formulated the gap equation for quartetting in 
fermion systems in analogy to the BCS gap equation. 
The mass operator of quartet-BEC with effective four- 
body vertices and in-medium four-body Schrodingcr 
equation was derived with the Dyson equation approach 
to correlation functions. However, the full expression of 
the quartet mass operator is too complicated to be eval- 
uated numerically in a direct manner. The biggest prob- 
lem stems from the many dimensional integrals over mo- 
menta. We, therefore, introduced some reasonable ap- 
proximations reducing the complexity considerably. A 
feature of great help is that the final answer is indepen- 
dent of the strength of the vertex function, making our 
approximation probably quite reliable. In our calcula- 
tion we also applied the same mean field ansatz projected 
on zero total momentum which was already so success- 
ful in our previous calculation of the critical tempera- 
ture [2]| . This feature, of course, reduces the numerical 
effort tremendously, since only a single OS wave function 
has to be determined selfconsistcntly by iteration. 

In this pilot work with an application to nuclear 
physics, we showed results only at zero temperature, how- 
ever, the formalism we presented is at finite temperature. 

We think that the results are of general validity show- 
ing qualitatively very distinct features from the pairing 
case. For example no well defined quasiparticles occur in 
the case of quartets. This is due to the fact that the quar- 
tet order parameter in the single particle mass operator 
goes along with the level density of three uncorrelated 
holes. Only the total momentum of the three holes is 
well defined and equal to the time reversed momentum of 
the incoming particle. The relative momenta of the three 
holes have to be integrated over, yielding a strong imag- 
inary part of the mass operator smoothing out any indi- 
vidual single particle structure. Another remarkable fea- 
ture already encountered in our previous work in Ref. [2]| 
is that the self consistent solution seems to exist only 
from negative chemical potential until around zero, i.e. 
from the BEC, or strong coupling region until crossover 
region. Once one goes to positive /i's the solution breaks 
down. This effect is even more pronounced here than it 
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was in our study [21| for the critical temperature. This 
can be traced back to the fact that the 3h level density 
goes through zero at 3/i for > 0, that is just at the 
place where the quartet correlations should occur. Actu- 
ally this feature is present for all multi-particle multi-hole 
level densities. The only exception being the single par- 
ticle level density which is finite at the Fermi level. This 
unique feature makes that pairing also is unique and for 
instance allows for a weak coupling situation with a co- 
herence length of the pair orders of magnitude larger than 
the interparticlc distance. It remains an open problem, 
how to formulate a more general theory which continu- 
ously goes from the quartet case to the pairing case. 
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Appendix A: Dyson equations for 
multi-particle-hole Green functions 

We shall shortly review the basic formulation used in 
the present work. We here extend the approach for real- 
time Green's functions (GF's) at zero temperature in 
Ref. [3^ to real-time GF's at finite temperature (25l[3l|. 

The Hamiltonian in a fermion system with two-body 
interaction is 

K = H-fiN = T + V-fiN 

= ^£lc\ci + ^^Vi2,34c\clc4,C3. (Al) 
1 1234 

The ci, c\ are fermion annihilation and creation opera- 
tors with an arbitrary quantum number 1. 

A real-time multi-particle-hole GF at a finite temper- 
ature is defined by 

^ ( -i{T{A^{t)Alit'))) (chronological), 
\-ie{t-t'){[A^{t),Al{t')]±) (retarded) ^ ' 

where (...) means the thermal average, T is the time 
ordering operator, A„ is an arbitrary operator consist- 
ing out of individual fermion operators ci and c\, and 
is the anti-commutator or commutator. The 
time dependence of the operators is given in the Heisen- 
berg picture Aa{t) = e'^^A^e-^^*. 

Note that, although we treat chronological GF's below, 
the change to retarded, advanced and Matsubara GF's 
goes as usual psj . Below, when we go from time space 
to Fourier space, we always will go over to retarded GF's 
without mentioning it explicitly. 



The superscript (ipjh) in Eq. (jA2|) means i- 
particle j-hole GF, where i (j) is the number of 
the annihilation (creation) operators in A^, e.g. 

G^^f}it-t') = -^{Tic,it)c\,it'))), G'^X^],,{t-t') = 

-i{T{{clci)t{c\,C2')t')) with (4ci)t = 4(i)ci(t), etc.. 
The Dyson equation for the ipjh GF is [30] 

+ J2[ dt"M^,p {t - t")G^;^^ it" ^t'), (A3) 

with 

= 5]([K,y]_,4]±)AA-i, (A4) 

l3 

Ma;a'(i-i') 

= -z5](T(K(t),F]_[y,4(t')]-))in..A/'^-:^„(A5) 

A/'a;a' = (K,4']±) (A6) 

where £„ is defined by [Aa,T — ^N]^ = Eq^q, and A/""^, 
is the inverse of the matrix Afa-a'- 

According to the time dependence, we shall call M^.^, 
static mass operator and Ma-a' [t — t') dynamical mass 
operator. 

Since [A^ (t) , _ is the operator A^ augmented 
by one annihilation and one creation operator, 
(T( [A„ (t) , _ [V, Al [t")] _ ))i„. becomes a (i + l)p(j -f 1 )h 
GF in the dynamical mass of Eq. (jA5[) . The index 'irr.' in 
Eq. (jA5[) stands for the mass operator being irreducible 
with respect to a cut of zpjh lines (26j . 

From the Dyson equation of Eq. (jA3p we can see that 
the Fourier transform of the bare ipjh GF is 

G°j;r^H = ;^- (A7) 

Appendix B: Quartet mass operator 

Here we derive the quartet mass operator of Eq. ^ 
shown in Sec. IIIIl 

Notice that we shall use summation convention for 
repeated indices and neglect all terms except the ones 
which are associated with the quartet order parameters. 

From Eq. (|A3|) . the Dyson equation for the Ip GF 
is [2^ 

+ J dtM-f ih - i3)4;2 ^ (^3 -t2). (Bl) 
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The Fourier transform yields 

(a. - e,)G^l.^[io) = <5i2 + M\^^\u)G[f {u^) (B2) 

with 

Therefore 

= G%^\u) + G%^\u:)M[^!\u)G^l^{ul (B4) 

where 

5^o 



(B5) 



UJ — El 

For the Dyson equation of the 2plh GF one obtains 



M 



(2plh) 



, (u;)G^^^]^'^ , (w), (B6) 



a^a^,Z3;a2a^,Z2 



where 



1 



(3p2h)irr. 



h\h\h'l ,yxv\-Mh'^h'^ .yiy'^ 



1 



a3a',Z3;a2a',Z2 



and 



-''^aia'j ,2:i;a2a^,2:2 ~ ( [''zi '^a'j^ , cj^^ ''a^ ''^2 ] + ) • 

Therefore 



(B7) 



(B8) 



^(2plh) 



^ ^0(2plh) 



ai ,zi ;a2a2 ,Z2 

' (^) 

ai a-j^ ,zi ;a3a2 ,2:3 ^ ' asa.^,z:i;a4a^,Z4 

xM^'P,''^' , (w)G^'P^'^^ , {u;), (B9) 

a4a4,Z4;a5a^,Z5 \ ' 050'^ ,25 ;a2a; ,Z2 ^ /J \ / 



^5ti^,Zi;a2a^,Z2 ^ 



with 



0(2plh) (^)^ A/;.a;,z,;a2a^,.2 



and 




FIG. 9: Graphical representation of P'^^ in Eq. (|BT3)| . The 
dots is represent the two-body interaction. 



where eij,,, = e 



+ •■• in (IbToI) . 

Substituting the 2plh GF of Eq. ()B9[) into the mass 
operator of Eq. (|B3p we obtain 



r(3) 



(3p2h)irr. 



^yiy'i'Mb'^b^^ ,y2y'2 



where 



xF 



.(3) 



(3)* 

2y2y2-b2b'^b'. 



(B12) 



lyiy[,bib[b'^ 



-Vlz^a^a-G^^^^r' , {u:)N-\ 
n -Lzi,uiu^ aia-^,zi\a3a^.Z3^ ' asa^ ,Z3 ;a4a^ , 



2^aiyi,bib'-^^ziy[5a'^b" + 2 ^i4y'i '^^li'i '^i4&" 



' r)'^Zib'{,yiy[^aibi5a'J,'_^ 



(B13) 



In (|B12p we omitted the term derived from the first term 
at r.h.s. of (jB9p because it is disconnected with the quar- 
tet order parameter. In the effective three body vertex 
r(3) of Eq. (jBTSl) we give not the exact 2plh GF but the 
free one, as we only consider in this work the lowest or- 
der approximation, though the exact 2plh GF figures, in 
principle, in the right vertex F*^"^^ of Eq. (|B13[) substitut- 
i ng di ll) into In Eq. ((BT3| . using Eqs. ((BT0| and 

(jBlip . we obtain 



a\a-^.z\:a:^a^,z^ ^ ^ c 
^aia^ ^a'-^a'^ ^ziz^ 
^ ~ a' H- ^z^ 



CL^(l^,Z^;CL4CL. ,Z4 



(B14) 



and thus Eq. (|B13|) is consistent with Eq. (fTTj) . 
Furthermore, the Dyson equation for 3p2h GF is 



bl b[ b'-[,yiy[ \b2b'^b'^ ,^21/2 



(3p21i) 

61 b'l b" , yi y 1 : 63 63 63 , t/3 ^3 



xG 



(3p2h) 



b^b'^b'^ ..ysy^Mb'r^b'^ ,y2y'2 ^ 

The mass operator is then given by 



(B15) 



M, 



(3p2h) 

bib[b'^ ,yiy[;b2b'^b'^ ,y2y'2 
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2 "^bix'-l ,cic'{'^yixi 5y'_^x'-^ '^fc'ic'j ^b'-[c'l 



~^ T^h'lx'l ,c'lc'l' "Vixi Oy'^x[ "hi ci Ofc'^ c[ 
~ '2^yic'l' ,xix'-[ ^y'^x'-^^bici ^b'-^ c[ ^b'^ c'-[ 
~'2^y{c'l',x'-^^x"^yixi5bici^b'^c'-^^^b'-[c" 



xG 



i(4p3h)irr. 



"^b^x'^ ,C2c'^' ^y3X2^y'^x'2^b'^c'.^^b'^c'^ 



2 
1 

1 



+ ^Vbir^'^ic'^c'.^' 5y^X2^y'^x'Jb3C25b'^c'^ 



yaC.^' ,X2x'^"y'^x:-^ 



Sy'x'Sb3C2Sb'c'Sb''c:;^ 



1 



'Vy'r^c'.^' ,x'^x'^5y3X2^b3C2 ^b'^c'^ '^''3 ^2 



and 



bs b!, b!,' ,y 3 ; 62 62 ''2' > ' 



-^fci ti'i f)" , J/i J/J ;fc2 bj, ''2' ^2/2^2 

J J J 



(B16) 



Therefore, one obtains for the 3p2h GF 



^(3p2h) / X 



G 



,0(3p2h) 

bi b'j b'l,yiy[;b2 b'^ b'^ ,1/22/2 
,0(3p2h) 



^bib'jb'/,j/ij/J;b3b^b!,',j/3j/^ ^^^'^b^b'.-^b'^ ,y3y'^-bib'^b'l ,yiy'^ 



XM 



(3p2h) 

bib'^b'^ ,yiy'^-br,b'^bi/ ,ysy'^ 



xG'(3p2h) , ^ 

b^b'^b'^ .y^y'^-Mb'^b'^ ,y2y'2^ 



(B18) 



where 



G 



,0(3p2h) 



bi b'j b'j' , y 1 y ^ ; b2 bjj b!j' , y2 ^2 
•M)i b'j b'/ ,yi y ^ ;b2 b^ b^' ,2/2^3 



(ebib'jb'Z 



(B19) 



'2/1^1- 



Substituting 3p2h GF of (jBlSp into the mass operator 
of Eq. (|BT2)) leads to 



= r 



'(^^ (w) 

,(4p 

CiC 



^g(4p3h)irr. ^ ^ ^^^^^ , , „ (w) 



xF 



(B20) 




FIG. 10: Graphical representation of P'*' in Eq. (|B2TJ. 



with 



F 



(3) 

iyiy[,bib[b'^ 



'b'.'H 



^g,0(3p2h) 

bib'jb'//yiyj;b3b^b^',y3y^ I'' b^b'^b'.^^ ,y3y'^-J>ib'^b'l ,yiy'^ 



'jVbix'^.ci c'/' ^y4Xi ^y'_^x[ ^b'^c[ ^b'lc'( 
+ 2 ^b^a;'/ c'/' <5y42;i ^y\x\ ^bici ^fc^'c'/ 
+ Vb'^x'l , c'{ c'l' ^yiXx Sy'^x[^b4Ci S),'^ c[ 



' 2 "a4c'/' dy'^x[ObiCi 5b'^c\ ^b^'c'/ 

' 9 ^aiC'i" ,an'ix'j' Sy^xi SbiCi Sb\c\ ^b'^c'^ 



(B21) 



Here we also omitted the terms associated with the first 
term at r.h.s. of Eq. (jBlSp because it is disconnected with 
the quartet order parameter. Besides we introduced the 
free 3p2h GF in F^^) of ([MT]) fo r the same reason as 
this was done in Eqs. (jB12p and (jBlSp . The graphical 
representation of F^"*^ is shown in Fig. 1101 

The 4p3h GF can, for our purpose of a-particle con- 
densation, approximately be decomposed into the order 
parameter and the free three body GF as follows (in anal- 
ogy to what is done in the pairing case) 

G^^Pf^i „(h-h) 

CiC-^C-^C^ ,XiX^X^ \ C2C[^C2C4 ,X2Xi^xi^ ^ -L / 

t t t \ / t t t t 



— ^^^^Xx^x'^ ^x'l "^c'/' Cc'/ Cc'^ Cci ) 1 1 (Cj,^ C^^ C^n C^n, Cx'^Cx'^Cx2)t2) 

x(?'(4iC^;C^,;')ti(c4'C4C^2)t2) 

= (Cc'/' Cc'/ C,; ) (c\^ c\,^ c\n c\,n ) 

xG°(^^) (tl-i2) 

X\X-^X-^\X2X\^X\^^ ^ ^> 

The Fourier transform is 



(B22) 



G 



(4p3h) 



12 



(cc;"Cc;'Cc;c,J(4cJ,cT„cT,„} 



t t t t 

'~* ^ c ^ ^ 

The free 3h GF is explicitly given by 



(B23) 



(o0(3h) , , 

Xix'^x'^ •^X2x'r^x'^ y / 

fxifx[fx" "I" fxifx[fx" 



uj + e 



with Pxix[x'{-X2x'^x '^ as in Eq. dH). 
Therefore, Eq. (|B23[) becomes 



G^^P?^', (o.) 

(cc;"Cc;'C,;Cci)(4c^,c^,,cJ,„) 

fxifx'^fx'^ ^~ fx^fx'^fx'-^ 



F^x X x'^x^^ ' x^x^^xi^ ' ( 25^ 



Finally we obtain for the quartet mass operator, 

Ixix^x^ ;ciC-^^C-^ \ / \ <--|^ c-j^ L-i / 



^ fxifx'^fx'^ ^" fxifx[fx'^ p 



><(44c^.ct )r^,^„^^^ (B26) 



Appendix C: in- medium quartet order parameter 
equation 

We here give a short derivation of the in-medium quar- 
tet order parameter equation. It is derived from the 
Dyson equation for 4p GF in static approximation of the 
mass operator. The Dyson equation for 4p GF derived 
from Eq. (|A3p without dynamical mass operator is 

~ ^1234^ G'^23i.i'2'3'4'(^ ~ 
= ■A/i234;l'2'3'4' 

4"Gl234;l'2'3' 

with in the ladder approximation 



+ ^1234:l"2"3"4"G^234:l'2'3'4'(* ~ (^^1) 



14234;1'2'3'4' 

([[C4C3C2C1, cJ„4„c35,,4„]_)7\/'~2./3"4'';l'2'3'4' 
1 



(1 - Pl - P2)l'l24'2"533'(544' 



+ 2^"'^ ^ Pi ^ P3)^'l3a'3"522"544' 
+ 2^-'^ ^ Pi ^ P4)l'l44'4"522"533' 

+ permutations 



(C2) 



where we approximated the correlation functions by fac- 
torizing them into products of single particle occupation 
numbers. 

The Fourier transform of 4p GF is 

^(4p) / ,N 

'-^1234;l'2'3'4'>>'^i 

~ ""1234;1'2'3'4' (.^i 

+ ^1234;1"2"3"4" ('^)-^l"2"3"4";5678 

X^5678;5'6'7'8'G5fg?7,g,.^, 2/3/4/ (w), (C3) 

A/i234:l'2'3'4' 



where 



r^0(4p) / X 

"'1234;1'2'3'4' \^ ) 



— £1234 



(C4) 



From the spectral representation of the 4p GF we shall 
only retain the ground state because of its condensate 
character. Therefore, at the ground state pole, i.e. at 
a; ~ 0, we obtain 

(C4C3C2C1} = ^V^1234:1'2'3'4'(C4'C3'C2'C1'}. (C5) 

£1234 

This equation corresponds to Eq. Q in the pairing 



Appendix D: Approximate mass operator in Eq. (|17p 

After the derivation of a single particle mass operator 
containing the quartet condensate in App. |b1 we easily 
recognize that its expression, for instance the vertex T^^'> 
in Eqs. ^ and ^ (or (|HT3)) and (|B2T|) ) it contains, 
is of considerable complexity, prohibitive for a direct nu- 
merical application, especially due to high dimensional 
integrals. We, therefore, will have to proceed to a quite 
intensive study of vertices, followed by reasonable ap- 
proximations, in order to reduce drastically the numeri- 
cal difficulty of the expressions. 

A first purely formal simplification which we will intro- 
duce, is to consider instead of the quartet case, only the 
trion case. Trions are fermions and one would have to 
develop a whole proper philosophy to introduce an order 
parameter for trions. However, we will ignore this diffi- 
culty here, proceed formally as if a trion order parameter 
existed in the same way as a quartet order parameter, 
and explain for this much simpler case our strategy. This 
can be done without loss of generality and in the end we 
simply will give our results for the quartet case which can 
be derived in complete analogy to the trion case. 

Let us, therefore, begin with the expression of the sin- 
gle particle mass operator containing a 'trion condensate' 
which analogously to Eq. (|B26p is given by 



Afl;l'(w) 

_ r(3) 



f'2i — ,/*3 
^123:456 (^) (C6C5C4) , . (^22' ^33' 



UJ + £23 



X(ci'4'4')ri'2'3';4'5'6'(w) 



(3)* 



<523"532') 

(Dl) 
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FIG. 11: The graphical expression of the first term on the 
r.h.s. of (fTTj). 




FIG. 12: Topologically same diagram as in Fig. 1111 



with r'^'^) the three body vertex already given in pT|) and 
derived in Eq. (jB13[) . The first and rather obvious ap- 
proximation we perform is to make r*-^-* 'instantaneous', 
that is cj-independent. A standard procedure for this is 
to put the vertex 'on the energy shell'. This procedure is 
not always defined unambiguously. In the present case, 
one possibility certainly is to put in F'^', a; = —£23, i.e. 
the energy where the mass operator (|Dip is resonant. 

Let us for the moment only investigate the first term 
on the r.h.s. of ([TT|) with this prescription. We obtain 



.(3,1) 
123;456 

1_ 

- ^^^13,3'6- 



-£23) 



-£2 



^3' — £6 



t'3'2,45 



(D2) 



A graphical interpretation of this term is given in Fig. 1111 
This graph also can be interpreted as a particular 
second order term of a three body scattering process 
(123)— )-(456), graphically represented in Fig. [T^ 
The intermediate propagator between the two vertices is 
given by (the energies are given by the prop- 

agators which are cut by the vertical line in Fig. [T^ . 
We, therefore, see that the static, on-shell part F(^^(a; = 
—£23) of (jDl[) can be interpreted as a second order three 
body scattering taken at w = 0, i.e. in reality at three 
times the Fermi energy. This certainly is a reasonable 
reduction of the second order process to a static vertex 

r(3,i)_ 

Proceeding in the same way with the second term on 
r.h.s. of ([TT|) . we arrive at a second order process anal- 
ogous to the one of Fig. [T^] as shown in Fig. [T31 We see 
that Fig. [T51 corresponds to Fig. [T^] with particles 2 and 
3 permuted. 

Let us now consider the third term in (1111) . In analogy 
to processes of Figs. [T^] and [13] this corresponds to the 
graph of Fig. (TD 




FIG. 13: The graphical expression of the second term on the 
r.h.s. of (HIJ. 




FIG. 14: The graphical expression of the third term on the 
r.h.s. of dnj. 



If we want to interpret this third term also as a static 
second order three body process with an intermediate 



propagator 



at cj = 0, contrary to first and 



-£l — £2/— £6 

second terms of (|lip . we have to use as on shell prescrip- 
tion in (fTTj) w = £1 -I- £4 -I- £5 -I- £6. We, thus, suppose 
from now on that all vertices in F'-'^^ are obtained from a 
second order three body scattering process at a; = 0. 

After these preliminary considerations, let us now in- 
troduce in a phenomenological intuitive way a single par- 
ticle mass operator containing a trion order parameter in 
a way analogous to the pairing case 

^Imri23;l'2'3'(^ - ^')G^3,2'3'(-^'), 

(D3) 

ImTl23;l'2'3' (w) 

-27r5(w)Vi23;456'/'l56V'l'5'6'^4'5'6';l'2'3', (D4) 



with 



and 



^123:l'2'3' 

-z{^\2.y2'hz' + Ul3,l'3"^22' + W23,2'3'^ll')i (D5) 



where we abbreviated the the order parameter as 



(C3C2C1) = -0 



123- 



(D6) 



The 'three body T-matrix' ((D4)) in ((D3)) is contracted 
with the antisymmetrized two hole propagator 

G^3,2'3'(-^) = -^^^^^('522'fe' - <^23"532'). (D7) 
W + £23 



With the above definitions (jD3j)-(|D7| we indeed see 
that the 'trion' mass operator is constructed in full anal- 
ogy to the one of the pairing case in Eq. It is evident 
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that in the quartet case we would define a four-body T- 
matrix as 



ImTi234; 



1'2'3'4 



(w) = -27r(5(tj)Ai234At' 



with 



and 



Ai234;5678 = Vi234;5678 (C8C7C6C5) 



(D8) 



(D9) 



1^1 



1234;1'2'3'4' 



--(;i2;l'2"533"544' + ^'yi3;l'3'^22"544' 



-permutations 



(DIO) 



and contract it with a 3h GF. 

Let us now investigate whether we can make contact 
of the mass operator in (|D3p with the one of (|D1[) . To 
this purpose, let us use in (|D4p the equation for the trion 
order parameter which in analogy to (jC3p is given by 



V'123 = — 



1 

£123 
(1 



-(1 



(1 

Pi - 

P2 - 



Pi 



P2)-jVi2.V2'h:i' 



Ps) 2^13,1'3"522' 
P3) 2^23,2'3'l5ll' 



V'J'2'3' (Dll) 



and investigate for the moment only the expression cor- 
responding to the second term on the r.h.s. of (jPSp . We 
obtain 



1 



2«13,46^25^l56 



2''^13,46'^25- 



a23 
1 



(1 - P4 - Pb)-jVi^A,w5&& 



+ {l - Pi- Pe)-VA(,A'&555' 



P6)2^56,5'6"^44' 



V'4'5'6' 



(D12) 



We notice some similarity with the second order vertex 
(|Dip . The difference stems from the presence of the oc- 
cupation numbers pi in (|D12[) and from the fact that 
there is more than one term. In principle the occupation 
numbers are the correlated ones. In the main part of the 
paper, we have seen that our theory practically only is 
valid for /i < implying that the occupation numbers 
remain small, see Fig. [51 We, therefore, can neglect the 
occupation numbers in (jD12[) to good approximation (for 
the p, < regime). Then, we can read off a vertex from 
(|D12p of the following form 



r 
1 



:3,2) 

123;456 



1 



7^13, 3'6 

t — £23'6 

1_ 1 

h-Wl3.46' 

4 -£246' 



1_ ^ 

^'3'2,45 + T''^13,3'6'025- 



-£23'6' 



-V3'6',46 



-■^26', 56 



(D13) 




FIG. 15: The graphical expression of the third term on the 
r.h.s. of (|DT3)| . 



inproper term for a vertex (see discussion below). The 
third term can be graphically represented as shown in 
Fig. [TS] We see that the term of Fig. [15] is obtained from 
Fig.[T2]in permuting indices 4 and 6. It is thus equivalent 
to the first term in (jD13[) . This game can be repeated 
in the same way for the first and third term in (|D5p . 
The conclusion is always the same: there are twice as 
many terms as equivalent terms from the on shell version 
of (jlip and in addition each time, there appears a dis- 
connected term which should not be present in a vertex. 
This latter problem is, however, not a real one. The order 
parameter is fully correlated and no disconnected pieces 
can appear. The term which is disconnected in (jD12[) 
will certainly become connected in higher orders. There- 
fore such disconnected terms only serve to renormalize 
the vertices of the connected terms. In the end all bare 
vertices in the connected diagrams should be replaced by 
two body T-matrices and no disconnected terms would 
appear anymore. We will not further dwell on this point 
and simply discard the disconnected terms. 

The conclusion, therefore, is that from (|D3p we get 
twice as many terms for the vertex F'-'^-' as is given in 
Eq. (|TT|) . Otherwise the terms are equivalent under the 
condition that in pT|) we take the above discussed on 
shell prescription. The correct procedure, therefore, is to 
divide the vertices Vi23:i'2'3' in (|D4p by a factor two. In 
the quartet case we have to divide the vertex in (IDIOI) by 
a factor of 8. 

In our numerical application (quartet case), instead 
of calculating all the 36 terms, where several ones give 
equal contributions, resulting from squaring the vertex 
(jPlOP where some are more difficult to calculate than 
others, we take a very pragmatic point of view and keep 
only the first term on the r.h.s. of (|D10p simulating all 
the others by a factor A' as in Eq. ([T5|) . In doing so, 
we suppose that all terms have more or less the same 
analytic structure. This is certainly the case, since all 
the terms are dominated by the behavior of the 3h level 
density whose typical structure is displayed in Sec. [V] It 
also is fortunate that the final result does not depend on 
A', as we now will demonstrate on the example of pairing. 



Appendix E: Independence of the parameter A' in 
Eq. ([T8|) 



It is realized that the first term is equal to expression 
(jPip . The second term is disconnected and, thus, is an 



All arguments and derivations for the three-body ver- 
tices can directly be generalized to the four body case. 
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In view of the fact that our favorite expression ([8]) is too 
complex for a direct numerical realization, we take a very 
pragmatic point of view for a pilot application. We have 
seen that some terms appear in all three forms we dis- 
cussed for the vertices. We shall take one of those terms 
for numerical calculations with a parameter A' in front 
of the vertex which should mock up the influence of fac- 
tors and also of topologically different graphs. However, 
we want to stress again that we do not think that topo- 
logically different vertices will finally give rise to different 
analytic structure of the mass operator. Essentially these 
additional terms will again only renormalize the vertices. 
The form which we then retain is just the first one on 
the r.h.s. of (DIO), since it is the one where the single 
particle motion directly couples via an interaction to the 
quartet amplitude, that is 



P 



uj- E' oj + E' 



{cc) 



1 - 2p' 

2£~ 



v(cc) 



(E3) 



with E'^ = £2 + D'^^ D' = v'{cc), and p' = {1 - s/E')/2. 
The equation for order parameter then reads 



(E4) 



where we should pay attention that the v in that equa- 
tion is still the original one, as in our four-body order 
parameter equation of (|C2[) and (jCSp . 

From above we obtain the gap equation for A' 



Ai234 A'-Wi2,l'2'(533'<544'(ci'C2'C3'C4') 



(El) 



With the separable form of the two body interaction [2l| , 
we then obtain expression ()21|) . As mentioned already, 
we will approximately absorb factors and all the other 
terms in renormalizing the vertex by a constant factor A'. 
This shall be a quite valid procedure, since the analytic 
structure of all the other terms is very similar. Fortu- 
nately our results will not depend on the strength of the 
rcnormalization factor A' of the vertex. This statement 
may be surprising but can be explained at hand of the 
standard BCS example as follows. 

Schematically, we write expression for single particle 
GF in BCS approximation of Eqs. ([T]) and ([2]) as follows 
(with self evident notation) 



A'^-E-^ (E5) 
Multiplying this equation with A , we see that E' ^ E 
and; therefore, nothing is changed in multiplying in BCS 
single particle mass operator the vertex by an arbitrary 
factor. Translated to our quartet problem, this implies 
that the question of multiplicity of diagrams discussed 
above is of no consequence for our numerical results. Of 
course, the addition of topologically different contribu- 
tions to the vertex F^^'' can change things slightly but 
again one can imagine that to good approximation it 
only renormalizes the vertex corresponding to the pro- 
cess which we treat explicitly. Therefore, our approxi- 
mate treatment of the single particle mass operator with 
quartet condensate should be quite safe. With our sepa- 
rable ansatz for the two body interaction, we then arrive 
to expression (jB20p of the mass operator. 



qBCS ^ 



1 



v{cc){c' c')v 



(E2) 

Appendix F: Preparation of mass operator for 
numerical application 



Now let us change in above equation vertex v to v' = X'v. 
Let us apply to all quantities containing this new vertex 
a 'prime'. Then BCS single particle GF is given by 



In this section, we describe the final expression for nu- 
merical calculation of the mass operator of Eq. (PT|) . 
We again give the approximated mass operator 



^quartet (^^^^^ 

d^k2 d^ks d^ki d^kv d^k2' d^kin d^k2" 
(27r)3 (27r)3 (27r)3 (27r)3 (27r)3 {2tt)^ (27r)3 

^ g-2(fei-fe.)V(4fc^)g-(fci.-£.0V(4fc^)e-(£i"-£2"')V(46^)(2^)3^(fc^ + fc^ _ fc^, 

X (^(|fci-|)<^(|fc2'|)'^(l^3|)'p(|fc4|)(2^)3(5(fci- + k2' + h + h) 
X (^(|fci-,|)(^(|fe"|)'^(|4|)<^(|^4|)(2^)3<5(A?i-, + fc2" + ^3 + h) 
^ f{k2)Kk3)f{ki) + f{k2)f{k^)f{ki) 

d^ko d^K d^k d^K' d^k' d^K" d^k" 



(27r)3 (27r)3 (27r)3 (27r)3 (27r)3 (27r)3 (27r)3 
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X ^{\^ + fc'D^d^ - fc'l + fc|)^(|| - fc|)(27r)35(A- + i?) 



X "^(l 



K" 



k"\M\ 



K" 



k"\)^{\ 



K 



/(fc2)/(|f + fc|)/(|f - fc|) + /(A:2)/(|f + fc|)/(lf 



fc|)(27r)35(i?" + A') 



(2^(2^' 



-2(fei + f 



+ i?|)/(|f + fc|)/(|f - fc|) + /(|fci + A-|)/(|f + fc|)/(|f - fc| 

, , , fe? , fcl-i? , 3if2 fc2 



(27r)3 



(Fl) 



The imaginary part of the above mass operator 
-ImM<i"^^"='(fci,w + i77) is 



-ImAf'5""'°*(fci,tj + i?/) 



d^if (Pk 
(27r)3 (27r)3' 



¥'(ly + fc|MlY-fc|)_ 



/(|fci+Al)/(| 



a: 



fc|)/(l 



K 



k\) 



+/(|^i+A''|)/(|| + fc|)/(||-fc|) 

ki-K 3K^ e ^ . 

Zvi m 4m m 



d^k' 
(27r)3 



XI 



A' 



A' 

y 



(F2) 



1. The fci = case 

The imaginary part of the mass operator for /ci = is 
jiven by 



(fK d^k 



(27r)3 (27r)3 



K ^ K ^ i2 
^(ly + fc|MlY-fc|)_ 

+ fc|)/(lf 



-/(AO/(lf + ^l)/(ly-fc|) 
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/(A')/(l|+Pl)/(lf 
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(F3) 
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Here, in the last equation, 



Therefore 



3m/i — muj 



maxO 



t = {p-K)/{pK), 
t' = [k' ■K)l{k'K). 



(F4) 

(F5) 

(F6) 
(F7) 



^2 

A:^ > 3m// - mw - y ^ kiK > (FIO) 

and 

< fc^ < 3m^ -rauj- — + fciX. (Fll) 

From the above equation, we obtain the hmits of integra- 
tion with respect to k: 



2. The fci / case 

On the other hand, -ImM'i^'^''*^* for fci 7^ is repre- 
sented by 



-IniM^"^'*<=*(fci,w-t-i77) 

1 pOO p 1 

dRK"^ I dti / dkk^ / dt 
1 Jo ./-I 



(27r)4 

^ g-2(fe:f + /fV4+feliftl)/b' 
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(27r)4 



1 nOO n 1 
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1 Jo J-i 



X e 



-2(fe^ + /f^/4+fei/fti)/b^ 



+/(|fci + i?|)/(|| + fc|)/(||-fc|) 
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kiK 
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{2n) 



kiK 

dk'k'^ I dt'e-^'"''" 



X ^(|-| + fc'IMI-f -fc'l)]'. 

Here, the fohowing condition has to be satisfied: 



1 < 



3m^ - mw - - - fc2 



< 1. 



(F8) 



(F9) 



/j2 3J^2 

Pmin{K) = \j max(3m^ - mcj - y kiK, 0) 



< k < y Smfi — muj — — 



Besides, from Eq. (jFlip . 



kl 'SK^ 



kiK ^ Pmax{K) 

(F12) 



fc2 

3m^ -~muj- ^ — + kiK > (F13) 



has to be satisfied. Therefore 

Pmin ~ max 



(ki — y' 9m/i — imuj , 



< AT < ^ I fci + \j 9TO/i - 3miLj - -y j = P„ 



3. Summary of integrals in ImAf'i"""=* 



(F14) 



With the Appendices |FT] and WM 

we can express 

TmA/'5""'°*(fci, w -I- ^77) by: 



If 9m/i — 3ma; < 0, 



(F15) 



Else: 

If ki = 0, 



-ImA/'i"'*'^*^*(fci,a;-fM7) 



2'!Tm 
2(27r)4 



dKK pe 



,.K ,.K 



-A'V(4fc") 
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, , d/fc'fc'^e-'^-"/^' / dt' 
(271-)^ Jo J-\ 



•p = \ 3mii — mu! 



3K 



'2 



PrnaxO = y -^i^im^L - muj), 

t = {k-K)/{kK), t' = [k' ■ A')/(fc'X)(F16) 
while if ki ^ 0, 

-ImA'/i"^'"='(A:i,w + M7) 
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-2(3m/j-mu;+fcJ/2-A:^/2-fc^)/b^ 



J{^j&m^, ^2mu-^- 2fc2)/(|| + fc|)/(|| - k\ 



K'^ K ^ K ^ 

+f{\l6mf, - 2mLU - — - 2fc2)/(|_ + k\)fi\- - fc|) 



1 



(27r)- JO 



dk'k'^-'^'"/"" 



dt' 



X '/'(l-f + fc'IMI-f -fc'l)]', 



Pmin{K) = ^/ max[3TO/i - - y kiK, 0], 
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p^g^^quartet calculated by Eq. 
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